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Abstract 

■ Nuclear many-body theory is used to study nuclear matter and finite nuclei at extreme isospin. 

In-medium interactions in asymmetric nuclear matter are obtained from (Dirac-) Brueckner theory. 
Neutron skin formation in Ni and Sn isotopes is investigated by relativistic mean-field calculations 
in DDRH theory with density dependent meson-nucleon vertices. Applications to light nuclei 
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are discussed with special emphasis on pairing and core polarization in weakly bound nuclei. Ap- 
proaches accounting for continuum coupling in dripline pairing and core polarization are presented. 
Calculations for the halo nuclei 8 B, n Be and 19 C show that shell structures are dissolving when 
the driplines are approached. Relativistic breakup data are well described by eikonal calculations. 
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oo '. 1 Introduction 
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^ ■ Tremendous progress has been made over the last decade in extending our knowledge into hitherto 
unknown regions of the nuclear chart. The experimental achievements [|I| were accompanied by com- 
plementary developments in nuclear theory. Nuclear structure physics has now access to nuclei with a 
large variety of proton-to-neutron rations, hence allowing to study nuclear forces at extreme isospin. An 
equally important new aspect is the strong reduction in separation energies close to the driplines. Under 
such conditions the conventional binding mechanisms known from /3 stable nuclei are likely to cease 
Ph . to be valid. There are experimental indications - which are strongly supported by theory - that single 
particle shell structures become dissolved in the vicinity of proton and neutron driplines. If confirmed 
this observation will have important consequences because it means that the seemingly well-established 
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$h ' prevalence of static mean-field dynamics is being replaced by dynamical processes from nucleon-nucleon 
— ■ (NN) interactions. In a very clean way this type of dynamics is observed in halo nuclei while in medium 
and heavy nuclei the dominance of the nuclear mean-field seems to survive. But with increasing neutron 
excess the character of the mean-field is dramatically changed by the strong enhancement of isovector 
interactions. For very neutron-rich nuclei theory predicts the formation of rather thick layers of almost 
pure neutron matter representing a genuine state of de-mixed proton and neutron fluids. 

Descriptions of nuclear interactions and nuclei far off stability necessarily afford strong extrapola- 
tions. Since data are still scarce the theoretical calculations rely on nuclear models, to be justified (or 
falsified) by experiment. An approach attempting a self-contained description of stable and unstable 
nuclei by nuclear many-body theory is presented in the following sections. 



2 Interactions in Strongly Asymmetric Matter 

An appropriate approach to in-medium nuclear interactions is Brueckner theory. Especially the relativis- 
tic effects accounted for with Dirac-Brueckner (DB) calculations have proven to describe nuclear matter 
properties rather satisfactorily. In ref . DB Hartree-Fock (DBHF) theory was used to investigate inter- 
actions in infinite matter over a large range of asymmetries. The Groningen free-space NN-potential was 
applied. Effective in-medium meson-nucleon coupling constants were extracted for the isoscalar a and 
tjj and the isovector p and S mesons, respectively 0. In all meson channels a pronounced dependence 
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Figure 1: Equation of state for symmetric nuclear matter and neutron matter obtained with 
the D3Y interaction (full lines), the relativistic DDRH interaction Q (dashed-dotted lines) 
and from variational calculations || (full squares and circles). 



Hence, to a very good approximation in-medium strong interactions remain intrinsically independent 
of isospin thus conserving a fundamental symmetry. Using the density dependent relativistic hadron 
(DDRH) field theory [[| the DB coupling constants have been applied to finite nuclei [|J, hypernuclei 
|| and neutron stars ||. DDRH theory is a microscopic formulation of relativistic mean-field theory 
J7| but describing in-medium effects by density dependent coupling constants rather than by non-linear 
meson self-interactions. The DDRH calculations describe binding energies, separation energies and 
other ground state properties of stable and unstable nuclei rather accurately on a level of a few percent 
j|] showing that the gross properties of nuclear matter and finite nuclei are described very satisfactorily 
with microscopic interactions. 

A clear indication for the incompleteness of a ladder-type interaction model is found in a re-analysis 
of the DDRH results with a mass-formula showing that especially surface energy contributions are un- 
derestimated [|J. An important class of interactions not accounted for by DB theory are contributions 
from ring-diagrams which will become important at low densities. Also missing are three-body interac- 



tions ||. In order to account for the missing contributions we follow the approach of |T0[ and introduce 
a semi-microscopic effective interaction by multiplying the Brueckner results by an additional density 
dependent vertex functional where parameters are fixed in symmetric matter. Results for the equation 
of state of symmetric and pure neutron matter obtained with DDRH theory are displayed in Fig.|l] and 
compared to the non-relativistic density dependent D3Y parameterization |TU . 



An interesting effect, predicted by many mean- field models, is the appearance of neutron- and proton- 
skins when the driplines are approached. In Fig.|2] DDRH proton and neutron density distributions for 
the isotopic chains of 48_82 Ni and 10 °- 140 Sn nuclides are displayed [|J. At the proton dripline ( 48 Ni, 
100 Sn ) proton skins are found being mainly caused by the Coulomb repulsion. With increasing neutron 
excess neutron skins starts to develop which, for example in the Ni-case, reach a thickness of about 
A = 0.8 fm at the neutron dripline, obtained here for 92 Ni. In skin nuclei a core of normal composition 
is coated by a layer of almost pure neutron matter. Investigation of this very particular state of nuclear 
matter by transfer reactions at REX-ISOLDE energies is discussed in B. 

In most of the following calculations for dripline nuclei non-relativistic HFB and RPA (or QRPA) 
theory will be used. The main reason for doing so is that at present the extension of DDRH theory 
(as is also true for phenomenological RMF models) to a full dynamical description is still in progress. 
Relativistic interactions are well tested at the level of HF and HFB theory, i.e. in the limit of a static 
ground state calculations in the no-sea approximation. But only few is known about extensions to 
dynamical processes, which in a relativistic approach should take into account also the coupling to 
vacuum modes. 

The non-relativistic D3Y interaction is based on a parameterization of the G-Matrix in terms of 




Figure 2: Charge densities (upper panel) and neutron ground state density distributions 
(lower panel) for Ni- and Sn-isotopes The formation of neutron skins in the neutron-rich 
nuclides is clearly visible. 



relativistic Brueckner theory to reproduce adequately the empirical nuclear matter EoS an additional 
density dependent re-scaling was introduced in the particle-hole (ph) channel, enforcing agreement with 
the variational EoS of ref.|| of for symmetric nuclear matter. Using Landau-Migdal theory [I2| the 



appropriate (density dependent) ph interactions in the various spin and isospin channels are derived 
self-consistently. 



3 Pairing at the Particle Threshold 

The existence of 11 Li as a particle-stable system relies almost completely on the mutual interactions 
among the last two valence neutrons. Their low separation energy, S2 n =320 keV, points to the impor- 
tance of the mixing with unbound states. For a detailed description of the valence wave function the 
conventional BCS and HFB methods using representations in terms of mean-field wave functions are 
not suitable. Theoretically, continuum effects are properly described by solving the Gorkov-equations 
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(h-e+)&+ - A$_ = , (/i-e_)$_ + A$+ = (1) 

as coupled equations for the hole and particle components [|nj] which are denoted by $± and are coupled 
by the pairing field A. The single particle energies are e± = \ ±E where A and E (> 0) are the chemical 
potential and the quasiparticle energy, respectively. Note, that e± will in general not coincide with the 
eigenvalues of the mean-field Hamiltonian h. In this sense, the states Q± are off the (mean-field) energy 
shell and deviate from usual quasiparticle picture. 

A necessary condition for particle-stability is A < 0. Then, irrespective of the value of E, the 
hole-type solutions $_ are exponentially decaying for r — > oo. For E< |A|, an exponential asymptotic 
behaviour is also found for the particle-type components $+ and a discrete subset of eigenvalues E is 
obtained. For E> |A|, eq.(^) has to be solved with continuum wave boundary conditions for $ + leading 
to single particle spectral functions being distributed continuously in energy. Hence the quasi-particle 
picture, underlying BCS theory and, to some degree, also discretized HFB theory, is extended to a fully 
dynamical continuum description. 

Here, we consider only like-particle (S — 0, T — 1) pairing. Writing for protons (q = p) and neutrons 
(q = n) $+(r, e+) = u q (e + )F q (r,e + ) and <3>_(r, e_) = v q {e^)G q {r; e_), respectively, where the reduced 
wave functions F q and G q are normalized to unity and u q , v q correspond to BCS amplitudes, the pairing 
fields are then defined in terms of the anomal or pairing density matrices for protons and neutrons, 
respectively, 
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Figure 3: Singlet-Even ( 1 S , , full line) and Triplet-Even ( 3 Si, dashed line) S-wave phase 
shifts obtained with the D3Y-interaction in free space are compared to calculations with the 
Bonn-B NN potential (circles and triangles) |T6| . 



and the proton and neutron pairing fields are 

A 9 (r!,r 2 ) = Vss(ri,r 2 )Kg(ri,r 2 ) . (3) 

If A < 0, K q (r) and A(r) are guaranteed to decay exponentially for r — > oo. The pairing interaction in 
the particle — particle channel is described by Vse- It is still an open question whether the free space 
or an in- medium singlet-even interaction should be used (see e.g. 

Here, for the numerical calculations a local momentum approximation is introduced by averaging 
in momentum space over the non-locality and replacing Vse^i,^) — * MsE(kF(r))8(ri — r 2}- This 
corresponds to including the off-shell correlations into the strength factor M$e which, as a consequence, 
acquires an intrinsic density dependence. It is given by the in-medium on-shell SE scattering amplitude 
at the local on-shell momentum kp{r), i.e. the (two-body) S-wave component of the SE Brueckner G- 
matrix in the particle — particle channel. The density dependence is re-adjusted such that for p — ► 
the free space SE S-wave phase shifts are reproduced. A comparison to the Bonn-B phase shifts ([HJ is 
shown in Fig.[3[ 

In order to understand pairing in weakly bound nuclei in more detail it is instructive to extract 
from eq.(|l]) (non-local and energy- dependent) dynamical pairing self-energies in the particle and hole 
channel, respectively, 

E±(l,2;a;) = A(l)t<7 T (l,2;u;)A(2) , (4) 

where g±(l, 2; ui) is the mean-field single particle propagator in the complementary (hole or particle) 
channel. This leads to the equivalent set of (formally) decoupled integro-differential equations 
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Figure 4: Single particle spectral functions for s- and p-wave neutron states in n Li. The 
finite widths of the states is due to continuum coupling. The partial occupation numbers 
Nji are also shown, including d-wave contributions. 



showing that the generalized pairing approach may depart strongly from a conventional BCS description, 
especially when £± ~ 0{ e ±) becomes comparable in magnitude to the effective single particle energies 
e±. In the constant-gap BCS approximation A(r) — > — Ao the self-energies become 

£^ = -A ^ ; S BCS = -A ^ (6) 

Ujt Vji 

and eq.(||) reduces to the conventional BCS eigenvalue problem. 

Returning to the full problem, strength functions for 11 Li are displayed in Fig.|j. The strong deviation 
from a pure mean-field or BCS description is apparent by observing that besides the expected s- and 
p-wave components also ^5/2,3/2 strength is lowered into the bound state region. Remarkably, the mean- 
field does not support neither bound 2s nor 1^5/2,3/2 single particle levels and their appearance is solely 
due to pairing. The states acquire a finite width because of the coupling to the particle continuum. 

Theoretically, the proton (q = p) and neutron (q = n) densities in a systems like n Li are defined by 

P,W=E^/ A ^-<,(e_)|G 9J ,(r,e_)| 2 , (7) 
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from which the particle numbers 

N g = Y, (2J + 1) f de.v 2 qje (e.) = £ (2j + l)n qjl , (8) 



are found. The neutron partial wave occupation numbers Nji = (2j + l)rij£ are included in Fig.[|. In 
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Figure 5: Proton and neutron ground state density distributions (weighted by r 2 ) for n Li. 
The strong neutron halo component formed by s-,p- and d-states is clearly visible. 



are missing because of the smallness of |A|. In Fig.[|the proton and neutron ground state densities are 
shown, multiplied by r 2 in order to emphasize the differences in shape and the neutron halo component. 
Applications of the densities in high-energy elastic scattering of n Li on a proton target are discussed 
in Measured n Li response functions were analyzed by QRPA calculations in |]18|| . 

Rather than solving the full continuum pairing problem as discussed here a widely used approach 
is to represent the Gorkov spectrum by a set of discretized energy levels by imposing a box boundary 
condition on the solutions of eq.(|]) |L5[]. A comparison shows a good agreement between the full and 
the discretized calculations for binding energies, chemical potentials and ground state densities provided 
the box radius is chosen large enough (Rb ox > 50/m). 



4 Dynamical Core Polarization at the Dripline 

Dynamical core polarization is seen most clearly in nuclei with a single nucleon outside a core. Ap- 
proaching the driplines the core nucleus by itself is already far off stability containing weakly bound 
single particle orbits. Under these circumstances the surface tension and therefore the restoring forces 
against external perturbations are reduced. Such "soft core" systems are, for example, found in the 
neutron-rich even-mass carbon isotopes. A good indicator is the existence of low-energy 2 + states, 
decreasing in energy with increasing mass. Continuum QRPA calculations with a residual interaction 
derived in Landau-Fermi liquid theory from the D3Y in-medium interaction reproduce the systematics 
of 2 + states rather well. This leads to the conclusion that they are mainly of vibrational nature rather 
than due to static deformation as assumed e.g. in . 

The calculations predict a strong increase of the quadrupole polarizibilities in the carbon isotopes 
for increasing neutron excess with a maximum around 16 ' 18 C. This behaviour is due to the lowering of 
the first 2 + state which in 12 C is located at E x = 4.44 MeV and moves down to E x = 1.66 MeV in 18 C. 
The dipole (or eletric) polarizibility, however, changes only within 10% over the 10_22 C isotopic chain. 

In a system with core polarization the valence particle (or hole) obeys a non-static wave equation 

(H M F + Z pol (e)-e)y = (9) 

including the static (HFB) mean-field Hamiltonian Hmf and the non-local and energy-dependent po- 
larization self-energy S po ; describing the rescattering of the nucleon off the core thereby exciting it into 
states of various multipolarities and excitation energies followed by subsequent de-excitations back to 
the ground state (see Fig.^). 

During these processes the particle can be scattered virtually into high lying orbitals. The quantum 
numbers of the intermediate 2plh (or lp2h) configurations are only constraint by the requirement that 
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Figure 6: Diagramatic structure of mean-field (left) and core polarization interactions (right) 
of a nucleon in single particle state j. Interactions (meson exchange) are indicated by dashed 
lines. Core polarization leads to intermediate states {j'J c ) with a (QRPA) core excitation 
J c and a single particle state j'. 




Figure 7: Core polarization self-energies for low energy 19 C(l/2 + ) and 19 C(5/2 + ) states. 
Only the local parts (r = r 1 = r 2 ) are shown. 



respect to the inert core in its ground state. In nuclear matter and finite nuclei these processes are 
known to give rise to a depletion of momentum space ground state occupation probabilities |20], |2l| 
from the step function momentum distribution of a Fermi gas of quasiparticles interacting only by a 
static mean-field. In practice, the core excitations are calculated in QRPA theory thus extending the 
static HFB picture in a consistent way to a dynamical theory. 

In a finite nucleus the dynamical single particle self-energies S po ; affect separation energies and 
wave functions. In dripline nuclei the polarization self-energies are found to have a pronounced state- 
dependence. In Fig.[F] DCP self-energies for the lowest l/2 + and 5/2 + in 19 C are shown. In the 5/2 + 
channel the overall strength is small and a fluctuating shape is obtained. The l/2 + states, however, 
experience an additional attraction from an attractive surface- type self-energy which, in fact, provides 
the main source of binding for 19 C(l/2 + , g.s.). 

Theoretically, this amounts to use a multi-configuration ground state containing a single particle 
component - reminiscent of the static mean-field configuration - and a multitude of configurations where 
the valence particle is rescattered into other orbits by interactions with the core |TJ], 113. Coupling to 



the lowest 2 + and 3 states only as advocated e.g. in |26] cannot account for the complexity of the 
process. 

Results for energies and spectroscopic factors in the single neutron halo nuclei 11 Be and 19 C and 
the single proton halo nucleus 8 B are shown in Tab.[IJ. In n Be core polarization is causing the 



reversal of l/2 + and 1/2 states by supplying an additional attractive self-energy in the l/2 + channel. 
The single particle spectroscopic factor S n (l/2*, g.s.)=0.75 is in good agreement with recent transfer 



Table 1: Energies, spins and ground state spectroscopic factors from core polarization cal- 
culations. 



Nucleus 


r 


Energy keV 


S(f,g.s.) 


8 B 


3/2" 


130 


0.73 


n Be 


1/2+ 


510 


0.74 


17 C 


5/2+ 


760 


0.51 


S 19 C 


1/2+ 


263 


0.41 



Table 2: Theoretical one-nucleon removal cross section on a 12 C target and width of the 
calculated momentum distributions obtained with core-polarized wave functions and eikonal 
reaction calculations are compared to data [Q. One-proton and one- neutron removal reac- 
tions are indicated by (-lp) and (-In), respectively. 



Isotope 


Energy 


(7-1N 


FWHM 


exp 
C-liV 


FWHM exp 




MeV/nucleon 


mb 


MeV/c 


mb 


MeV/c 


8 B (-lp) 
10 B (-lp) 
17 C (-In) 
19 C (-In) 


1440 
1450 
904 
910 


104 
17.3 
124 
192 


75 
145 
132 
69 


98±6 
17±2 
129±22 
233±51 


91±5 
165±8 
143±5 
68±3 



are strongly suppressed as seen from the small spectroscopic factors in Tab .[I]. Hence, the valence 
nucleon is no longer attached to a definite mean-field orbital but exists in a wave packet-like state 
spread over a certain range of shell model states. In other words, mean-field dynamics have ceased to 
be the dominant source of binding. Rather, the 17 C and 19 C results indicate a new type of binding 
mechanism in dripline nuclei where shell structures are dissolved and binding is obtained from dynamical 
valence-core interactions. 

Core polarization also affects the low energy continuum region of dripline nuclei. The admixtures 
of core excitations may give rise to configurations where neither of the involved nucleons is in a state 
above threshold although the total energy is well above particle threshold. Such states are known as 
bound states embedded in the continuum (BSEC). They decay by coupling to energetically degenerate 
single particle continuum states. A known case is a 3/2+ state in 13 C which was predicted in [27] and 
observed experimentally |p8| . In 19 C BSEC-type configurations are indeed found in the s-,p- and d-wave 
channels below about 6 MeV above neutron threshold. They appear as resonances in the continuum 
strength functions. In Fig.|8| 19 C(l/2+) and 19 C(5/2+) DCP continuum spectral functions are shown. 
The BSEC structures are most prominent in the 1/2+ channel because the missing potential barrier 
will not support potential resonances. In other words, observing a 1/2+ resonance is a direct proof for a 
state with a non-trivial many-body structure. The 5/2+ state seen in Fig.|8| is part of a (much broader) 
conventional d-wave potential resonance being shifted into the low energy region by core interactions. 
Clearly, the existence of such close-to-threshold BSEC states will have important consequences for 
breakup reactions and, especially, for neutron capture in astrophysical processes. 

The theoretical overlap wave functions have been used in eikonal breakup calculations at relativis- 
ts energies. The FRS/GSI data for longitudinal momentum distributions |22], ^3| and removal cross 



sections [241 are well described as seen from Table 



5 Summary 

Nuclear many-body theory was applied to strongly asymmetric nuclear matter and the ground and 
excited states of dripline nuclei. Interactions were derived from Brueckner and Dirac-Brueckner calcu- 
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Figure 8: Continuum strength functions for 19 C(l/2 + ) (upper panel) and 19 C(5/2 + ) (lower 
panel). BSEC-like resonance structures are seen, especially in l/2 + channel. 



are not contained in the ladder approach. Pairing in weakly bound system was described in a approach 
taking into account the coupling to unbound continuum configurations. Dynamical core polarization in 
single nucleon halo systems was investigated by HFB and QRPA methods, also accounting for contin- 
uum effects. 

The overall features of dripline nuclei are rather well described. However, open and interesting 
questions remain on the structure of interactions especially in the low-density region. Here, the present 
treatment of interactions lacks full self-consistency by the semi-empirical adjustments to variational re- 
sults. Especially three-body interactions and induced 2-body interactions from ring diagrams in weakly 
bound asymmetric nuclear matter and finite nuclei have to be investigated further. Experimentally, 
this could be complemented by high-resolution measurements of spectral functions in the bound and 
the continuum region. 

The calculations strongly emphasize the importance of continuum coupling. They are found to 
dissolve shell structures close to the driplines and create BSEC structures in the low energy continuum. 
It will be interesting to investigate core polarization and continuum pairing for transfer and breakup 
reactions and capture reactions in astrophysical scenarios. 
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